Equilibrium solutions of the shallow water equations 
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A statistical method for calculating equilibrium solutions of the shallow water equations, a model of 
essentially 2-d fluid flow with a free surface, is described. The model contains a competing acoustic 
turbulent direct energy cascade, and a 2-d turbulent inverse energy cascade. It is shown, nonetheless 
that, just as in the corresponding theory of the inviscid Euler equation, the infinite number of 
conserved quantities constrain the flow sufficiently to produce nontrivial large-scale vortex structures 
which are solutions to a set of explicitly derived coupled nonlinear partial differential equations. 



The evolution of a fluid from a strongly random initial 
condition is generally characterized by one or more tur- 
bulent cascades of energy to larger and/or smaller scales. 
Whether energy flows to smaller scales via a direct cas- 
cade, or to larger scales via an inverse cascade, is deter- 
mined by a combination of conservation laws and phase 
space considerations. Generally, if only energy is con- 
served |l| (as for 3-d Navier-Stokes turbulence Q), its 
flow in phase space will be globally unconstrained and 
will spread out to arbitrarily high wavenumbers, even- 
tually draining all energy out of any large-scale macro- 
scopic flows initially present. Perhaps the most famil- 
iar example of this is the thermodynamic equilibration 
of a container of gas to a macroscopically featureless fi- 
nal state in which all energy eventually ends up as heat, 
i.e., microscopic molecular motion. If, however, one or 
more additional conservation laws are present (as for 2-d 
Navier-Stokes turbulence J 2J , or for deep water surface 
gravity wave turbulence p ] ) their multiple enforcement 
will generally not permit both conserved quantities to es- 
cape to small scales, and macroscopic structure, whose 
profile will be initial condition and boundary condition 
dependent, will survive. An example is the equilibration 
to a rigidly rotating final state of a gas in a cylindrical 
container with frictionless walls. The additional conser- 
vation of angular momentum along the axis of the cylin- 
der precludes a featureless final state. 

A long-standing problem has been the characteriza- 
tion of final states of systems with an infinite number of 
conservation laws [Q. Different values of the conserved 
quantities should then produce an infinite dimensional 
space of final states ||. The example of the equilibrat- 
ing gas motivates one to postulate that the macroscopic 
final state, be it featureless or not, should be thermody- 
namic in character, i.e., it should be an equilibrium state 
computable from the appropriate Hamiltonian using the 
formalism of statistical mechanics. In Q this approach 
was used to produce a full characterization of the equi- 
librium states of the 2-d incompressible Euler equation 
(the inviscid limit of the 2-d Navier-Stokes equation) , 
where the conserved quantities are the standard integrals 



of all powers of the vorticity. The equilibria were found 
to be characterized by a macroscopic steady state vor- 
ticity distribution ujq (r) obeying an explicit "mean field" 
partial differential equation whose input parameters were 
determined by the values of the conserved quantities. 

The 2-d Euler equation is the simplest of these systems 
in the sense that the incompressibility constraint V-v = 
reduces the dynamics to that of the single scalar vorticity 
field to = V x v = d x v y — d y v x , and the conservation laws 
then provide an infinite sequence of global constraints on 
its evolution. In this work we study a more complicated 
system of equations, the shallow water equations, an ex- 
tension of the 2-d Euler equation that includes a free 
surface with height field h(r) coupled to gravity g. The 
horizontal velocity v now becomes compressible (with 3- 
d incompressibility enforced via v z — -zV • v) but is 
assumed to be independent of the vertical coordinate z. 
The effective 2-d dynamical equations are: 



~Dt 



= <9 f v + (v • V)v = -gVh 
d t h + V • (hv) = 0, 



(1) 
(2) 



The first equation expresses the fact that the fluid ac- 
celerates in response to gradients in the surface height, 
and the second enforces mass conservation, i.e., the full 
3-d incompressibility. The Euler equation is recovered 
formally when g — > oo since height fluctuations are then 
suppressed. It is straightforward to verify that the ratio 
= ui/h is convectively conserved, Dfl/Dt = 0, imply- 
ing conservation of all integrals of the form 



(3) 



C f = / d 2 rh(v)m(v)} 



for any function /(s). These may be fully characterized 
by the function g(cr), — oo < a < oo, obtained from (|^) 
with f(s) — 8(ct — s), and g(o~)da therefore represents 
the 3-d volume on which a < fl < a + da. For general 
/ one then recovers Cf = / daf(a)g(a). Note that if 
u> = initially, then it must remain zero for all time. 
Initial conditions of this type then generate (nonlinear, 
in general) wave motions only f| . 
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The extension of the equilibrium theory to the shallow 
water equations is a significant advance because in addi- 
tion to the usual vortical motions they contain acoustic 
wave motions ||. The latter are known |hJ to have a 
direct cascade of wave energy to small scales. One then 
has the very interesting situation in which there are two 
competing energy cascades, and the question arises as to 
which one "wins." In particular, is it possible that the 
macroscopic vortex structures can "radiate" wave energy 
and disappear entirely? We will show that under reason- 
able physical assumptions a finite fraction of the energy 
remains in large scale vortex structures, and we will de- 
rive exact mean field equations for the equilibrium struc- 
ture. 

The statistical formalism proceeds in a sequence of well 
defined steps. First, the Hamiltonian corresponding to 
(111) and (0) is 



H= \j d 2 r(hv 2 + gh 2 ), 



(4) 



though the Poisson bracket yielding (|l|) and (H) from (^) 
is noncanonical [Q. Second, the partition function is de- 
fined as an integral over the phase space of fields h, v with 
an appropriate statistical measure. This so-called invari- 
ant measure is most easily computed if the dynamics can 
be expressed in terms of a set of variables, canonical vari- 
ables being an example, for which a Liouville theorem is 
satisfied. In this case invariant measures are any function 
of the conserved integrals, with different choices corre- 
sponding to different ensembles. In the Euler case || the 
field uj itself satisfies a Liouville theorem. In the shal- 
low water case no obvious combination of h, v or their 
derivatives meet this requirement. 

To circumvent this problem we transform to a La- 
grangian description, in terms of interacting infinitessi- 
mal parcels of fluid of equal 3-d volume, for which canon- 
ical variables are easy to construct. Thus, let a be a 2-d 
labeling of the system, and let r(a, t) be the position of 
the parcel of fluid such that, e.g., r(a, 0) = a. Since 
all parcels have equal mass, the conjugate momentum is 
p(a, t) — r(a, t) = v(r(a, t), t). The height field is simply 
the Jacobian of the transformation between r and a: 

h /h(r(a)) = dct(<9r/<9a) = d ai r 2 - d a2 r 1 , (5) 

where ho is the overall mean height. The Hamiltonian 
(0) now takes the form 



(G) 



while, 



w(a) = V x v = h(d a2 rid ai P\ - d ai rid a ,_Pi 

+ d a2 r 2 d ai p 2 - d ai r 2 d a2 p 2 ) 
q(a) ee V • v = h{d a2 r 2 d ai pi - d ai r 2 d a2 pi 

+ d a2 rid ai p 2 - d ai rid a2 p 2 ). (7) 



It is easily verified that the Lagrangian forms of ([j]) 
and (|J) follow from the Hamiltonian equations of motion 
f(a) = 8H/5p(a) and p(a) = -6H/6r(a). The Liou- 
ville theorem, which is a statement of incompressibility 
of flows in phase space, 

J2 I d 2 a[Sr a (a)/8r a (a)+Sp a (a)/5 Pa (a)} (8) 

= J d 2 a[S 2 H/Sr a {a)Sp a {a) - S 2 H/Sp a (a)Sr a (a)} = 0, 

then follows immediately and implies that the correct 
statistical measure is p{TL,{g{a)})Y\ a d 2 r{a)d 2 p{a). In 
the grand canonical ensemble, which we shall adopt, the 
function p is given by p = e~ /3lc , where (3 = 1/T is a 
hydrodynamic "temperature" and 



K = H 



dap(a)g(a) 



= U - J d 2 rh{v)p[Lo{v)/h{r)] 
= H-h { d 2 ap[u(a)/h(a)} 



(9) 



in which p(o~) is a chemical potential that couples to each 
level w(r)/ft,(r) = a. The partition function is now de- 
fined by 



Z[/3,M<r)}] 



an 



d z r(a) / d z p(a)e 



-0K 



(10) 



where N — * oo is the number of fluid parcels and AM is 
the usual classical delabeling factor. The thermodynamic 
averages of the conserved quantities are now obtained 
in the usual fashion as derivatives with respect to the 
chemical potentials, (g(o~)) = TSln(Z)/Sp(cr). 

One would now like to transform the integration in 
( |l0| ) back to physical Eulerian variables. The key ob- 
servation is that, from (0), f2 ee lo/Ii and Q = q/h are 
linear in p. Therefore, one may formally invert this re- 
lationship to obtain FJ a d 2 p(a) = Y\ a dQ(a)dQ(a)J[h], 
where, due to the particle relabeling symmetry (both 
V • v and V x v depend only on r and are then clearly 
invariant under any permutation of the labels a), the Ja- 
cobian J is a functional of the height field h(a) alone. 
The exact form of J will turn out to be unimportant. 
Similarly, (1/AH) J] a J d 2 r{a) = H a J dh{a)I[h], where 
I[h] is another Jacobian. The 1/Nl factor precisely re- 
moves the relabeling symmetry that, in particular, leaves 
the height field invariant. Finally, we replace the la- 
bel a by the actual position r, in which the equal vol- 
ume restriction on each fluid parcel implies that the in- 
finitesimal area of each parcel must be determined by 
dV = h(r)d 2 r — constant. Thus: 



Mll/^(a)/^(a) 



H / dh(r)J[h] / dfl(r) / dQ{r), (11 
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in which J[h] = I[h]J[h], and the mesh over which the 
label r runs is nonuniform and changes with each real- 
ization of the height field h. 

The statistical operator K, must also be expressed in 
terms of Q,£l,h. Only for the kinetic energy T — 
f cPrhv 2 does this require some nontrivial manipula- 
tions. Let the current j = kv be decomposed in the 
form j = V x ip — V0. One obtains then 



( hQ 



V x i Vx -V x i V 
V'^Vx -V-IV 



(12) 



The 2x2 matrix operator, which we shall denote Ch, 
appearing on the right hand side of (|lj) is self adjoint 
and positive definite, and therefore possesses an inverse, 
i.e., a 2 x 2 matrix Green function Qh(r,r') satisfying 
£hGh(r, r') = l<5(r — r'). An explicit form for Qh will not 
be needed. The kinetic energy is then T = | J d 2 rj ■ v = 
f d 2 rh{ipn + <f)Q), i.e., 



T = J d 2 rh{v) J d 2 r'h{r') 



(13) 



and the complete statistical operator is 

K = T + J d 2 rh{v) | l -gh{v) - M [fi(r)] j (14) 

The appearance of the factors h{v) and h(r') is crucial 
here because, as discussed above, dV — h(r)d 2 r and 
dV' = h{r')d 2 r' are both uniform for each given statisti- 
cal mesh. 

We finally come to the evaluation of the partition func- 
tion itself. This is accomplished with the use of the 
Kac- Hubbard- Stratanovich (KHS) transformation, which 
in discrete form reads for any positive definite matrix A, 



(15) 



where yi and Q may be vectors, and the normalization 
is M = y / det(27rA). This identity follows by completing 
the square on the right hand side and performing the re- 
maining Gaussian i nteg ral. We apply it to the discretized 
version of (0) and (|ll"| ) with finite dV, and the identifica- 
tions Aij = -p-^ruTi) @, Vi = pdVlSlixilQiixi)] 
and we introduce the notation (i = [^(xO, ^( x i)]- The 
continuum limit dV —> will be taken at the end. The 
partition function is now 



Z = ]J dh 



m 

l N[h] 



d^> l d<S> l / dQidSke?*, (16) 



where 



(17) 



in which is an appropriate discretization of the 

differential operator Ch- Notice that the inverse of Qh has 
led to the reappearance of the local differential operator 

At the expense, then of introducing the new fields ^, 
$ we have succeeded in producing a purely local action 
in which the integration over fij, Qi can be performed in- 
dependently for each i (for given fixed field h). However, 
we now arrive at a problem whose physical origin, as we 
shall see, lies precisely in the direct cascade of wave en- 
ergy. Thus, the chemical potential function /Lt(er) controls 
convergence of the Laplace transform-type integral 



JW[*i] = 



dQj 



(18) 



where (3 = f3dV corresponds to a rescaled hydrodynamic 
temperature T = TdV which is assumed to remain finite 
as dV — > — the object of this choice is to obtain the cor- 
rect control parameter for nontrivial hydrodynamic equi- 
libria in the continuum limit that, as we shall see, yields 
a nontrivial balance between energy and entropy contri- 
butions to the final free energy Hj. However, there is 
no corresponding chemical potential controlling Qi and 
the corresponding integral does not converge. Recall- 
ing that Q = (l//i)V • v, unboundedness of Q reflects 
unboundedness of small-scale gradients in the compres- 
sional part of v and in h fL2[ . Thus, taken literally, the 
direct cascade of wave energy leads to arbitrarily small 
scale fluctuations of the fluid surface that remain of fixed 
amplitude, i.e., a kind of foam of fixed thickness. Physi- 
cally, of course, such small scale motions are rapidly dissi- 
pated by processes that violate the approximations used 
to derive the shallow water equations, e.g., by some com- 
bination of viscosity and wave breaking |p_3|| . This leads 
to the following physically motivated assumption: dissi- 
pative processes that suppress wave motions lead to the 
interpretation J dQiexp((3Qi$i) — * <5(/3$j), i.e., to the 
vanishing of $i. 

With = 0, only the (1, 1) component of Ch con- 
tributes, and in the continuum limit dV — > the parti- 
tion function becomes 



where the Free energy functional is 
T= I d 2 r 



(19) 



(20) 



The key observation now is that (3 = f3/dV — > oo in the 
continuum limit. Thus, mean field theory becomes ex- 
act and equilibrium solutions are given by extrema of T . 
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This is why the integration over the field h(r), with its 
unknown Jacobian, is ultimately irrelevant. The under- 
lying assumption is only that the Jacobian is smooth, or 
at least less singular than e - ^, in the neighborhood of 
the extremum in the continuum limit. 

The extremum conditions ST/S^^x) = = SJ 7 /5h(x) 
yield then the mean field equations 



1 



Mr) 
[W(r 



2h(ry 



= h(r)W'[*(r)] (21) 
= W[V(r)]-gh(T), (22) 



By adding a source term J d 2 rh(r)T(r)Q(r) to JC, which 
serves only to replace \& by , 5 — r inside W, one may com- 
pute the equilibrium average (O(r)) = [5!F /8t{t)\ t =q = 
-/i _1 V ■ (/i _1 V*). It follows then that (j) = V X so 
that iS? is the stream function associated with the equi- 
librium current. Equation (|2l]) is in fact equivalent to 
il = — W'(^>), which guarantees that this is a true equi- 
librium solution satisfying O = 0, and equation (|2^) is 
equivalent to Bernoulli's theorem since it can be rewrit- 
ten as (l/2)v 2 +gh = 

As a simple example, in the case where = <7q 
over half the 3-d volume of the fluid and = on 
the other half, the chemical potential takes the form 



6(a) 



(To), and therefore by ( |is| ) 
Extensive numerical solu- 



2 M<0 

e 0W(s) _ e Pno _|_ e 0(ni-<7os) ^ 

tions for the Euler equilibria exist for this "two-level" 
system as a function of (3 and /Ui — /Uo- |]]7j]. In pre- 
liminary numerical work, we find that the shallow water 
equilibria generated by ( p2|) have very similar structure 
(with, for example, vorticity moving from the walls to- 
ward the center of the system as j3 decreases from positive 
to negative values), while the height field basically co- 
varies with the vorticity in order to maintain hydrostatic 
balance. Details of this work will be presented elsewhere. 

The techniques presented in this paper can be used to 
generate equilibrium equations for a number of other sys- 
tems with an infinite number of conserved integrals Q. 
The key insight presented here is that whenever such a 
system contains simultaneous direct and inverse energy 
cascades, the long time dynamics becomes very singular 
and additional physically motivated assumptions must 
be made in order to derive sensible equilibria. Our as- 
sumption, that dissipation acts to suppress the forward 
cascading degrees of freedom with negligible effect on the 
macroscopic state, presumably depends on the smooth- 
ness of the initial condition. Comparisons with detailed 
numerical simulations will be required to evaluate such 
effects. 

Note added: After completion of this work we be- 
came aware of an e-print [ fl4| where equations equivalent 
to (p2p are derived from a phenomenological maximum 
entropy theory. No statistical mechanical derivation is 
given, nor is the interaction between wave and vortical 
motions and the effects of waves on equilibration dis- 
cussed. 



[1] 

[2] 
[3] 
[4] 

[5] 

[6] 



[8] 



[9] 

[10] 
[11] 



[12] 



[13] 



[14] 



Real fluids are always viscous, but models of turbulence 
generally concern themselves with the "inertial range" 
where loss of energy due to viscous damping is small com- 
pared to that due to the cascade process, and an energy 
conserving model is appropriate. 

See, e.g., A. S. Monin, and A. M. Yaglom, Statistical 
Fluid Mechanics, Vol. 1, (MIT Press, Cambridge, 1971). 
See, e.g., G. Falkovich, V. L'vov, and V. E. Zakharov, 
Weak turbulence theory of waves. 

For several examples of such systems see, e.g., D. D. 
Holm, J. E. Marsden, T. Ratiu and A. Weistein, Phys. 
Rep. 123, 1 (1985). 

Study of these states is partly motivated by the con- 
straints they place on simulations of turbulent flow, e.g., 
the degree to which they properly preserve the conserva- 
tion laws. 

J. Miller, P. B. Weichman, and M. C. Cross, Phys. Rev. 
A 45, 2328 (1992). 

The statistical approach relies on the assumption of er- 
godicity of the dynamics. This assumption has been ex- 
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The linearized versions of these equations, dtv = —gVn, 
where n — h — ho with ho the mean surface height, and 
dtT) + ho V • v = 0, have longitudinal traveling wave solu- 



te 



(k-r-c|k|i) 



r)o(c/ho)ke 



i(k-r — c|k|t) 



with 



tions 77 

speed c = \fgh~o~. These waves become coupled through 
the nonlinear terms as the amplitude 7/0 increases. 
If u) = then v = — \7(j> leads to the pair of scalar equa- 
tions d t (f> = g{h - h ) + i|V0| 2 , d t h = V ■ (hV<j>). 
See, e.g., A. Balk, Phys. Lett. A 187, 302 (1994). 
Note that j3 < is required for positive definiteness of A, 
i.e., hydrodynamic equilibria often correspond to negative 
temperature states. This is explained in detail in ||. If 
(3 > one uses the the KHS transformation with Q — > iQ . 
Divergences in these gradients may in fact occur in finite 
time since the shallow water equations are believed to 
produce shock wave solutions. Appropriate continuation 
of the equations nevertheless allows the conservation of 
ui/h to be maintained even the presence of shocks (D. D. 
Holm, private communication). 

Dissipation processes also act on the vortex structures, 
but the conservation laws guarantee that microscopic 
fluctuations in Q, unlike Q, remain finite and hence will 
be dissipated much less strongly. 
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